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■ 1. Introduction 



ABSTRACT 



We use XMM-Newton data of the merging cluster Abell 3667 to analyze its metallicity distribution. A detailed abundance map of 
the central l.lxl.l Mpc region indicates that metals are inhomogeneously distributed in the cluster showing a non-uniform and very 
complex metal pattern. The highest peak in the map corresponds to a cold region, slightly offset South of the X-ray center. This 
could be interpreted as stripped gas due to a merger between a group moving from NW towards the SE and the main cluster. We 
note several clumps of high metallicity also in the opposite direction with respect to the X-ray peak. Furthermore we determined 
abundances for 5 elements (O, Si, S, Ar, Fe) in four different regions of the cluster. Comparisons between these observed abundances 
and theoretical supernovae yields allow to get constraints on the relative number of SN la and II contributing to the enrichment of the 
intra-cluster medium. To reproduce the observed abundances of the best determined elements (Fe, O and Si) in a region of 7' around 
the X-ray center, 65-80% of SN II are needed. The comparison between the metal map, a galaxy density map obtained using 550 
spectroscopically confirmed cluster members and our simulations suggest a recent merger between the main cluster and the group in 
theSE. 

Key words, galaxies: cluster: general - galaxies: abundances - X-ray: galaxies: clusters: individual (Abell 3667) 
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Cluster of galaxies are the largest virialized objects in the uni- 
verse. They form via gravitational instability from the initial 
perturbations in the matter density field. The cosmic baryons 
fall into the gravitational potential of the cluster dark matter 
halo formed in this way, while the collapse heat up the intra- 
cluster medium (ICM). At the high temperature measured in rich 
cluster, kT>3 keV, the ICM is highly ionised and its spectrum 
presents several emission lines, among which the most promi- 
nent is the Fe K- shell line at ~ 7 keV. As heavy elements are 
only produced in stars the processed material must have been 
ejected into the ICM by cluster galaxies. 

There is more and more observational evidence that various 
types of processes a re at work (e.g. ram pressur e stripping and 
galactic winds; see ISchindler & Diaferiol 120081 for a review), 
which remove interstellar medium (ISM) from the galaxies. 
Hence, the study of the metal distribution is a sensible way to 
better understand the thermodynamical properties of the diffuse 
gas and the past h i story of star formation in galaxy clusters 
(Arnaud et al . 1992; Renzini et al. 1993; Renzini 2004). 
Simulations have shown that the metal distribution in a clus- 
ter shows many stripes and blobs at the positions where the 
enrichment has taken place. Depending on the dynamical state 
of the cluster the metal distribution looks very different: in 
simple merger configuration (e.g. collision between two sub- 
clusters) pre-mergers have a metallicity gap between the sub- 
clust ers, post-mergers have a high metallicity between subclus- 
ters (Kapferer et a l]l2006b IKapferer et al.ll2007l ISchindlerl l2008: 
[Schindler & DiafericJ|200|)7^ 

Although to obtain metal maps from observation is not easy 
because a lot of photons have to be accumulated in each 



region to measure the metal abundance, several groups, us- 
ing XMM-Newton and Chandra data, have derived quite de- 
tailed metallicity maps (ISchmidt et al .120021: [Sanders et al.|2004|: 

" 0051: 

Werner et al.l 120061: ISanders & Fabianl 120061: lHavakawa etaD 
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Durret et al. 2005; O' Sullivan et al. 2005; Sauvaseot et al. 2005 
Wer ner et al.l 120061: ISanders & Fabianl 120061: lHavakawa et al. 
2006; Simionescu et al. 2009|). In general these 2D maps show 
that the distribution of metals in cluster is not spherically sym- 
metric, but it has several maxima and complex metal patterns. 
The range of metallicities measured in a cluster from minimum 
to maximum comprises easily a factor of two. 
Because clusters retain all the metals provided by their galaxies, 
X-ray measurements of the abundances of each element in the 
ICM enable us to examine the ratio of SN la and II. 
In this paper we present the results of the analysis of the clus- 
ter Abell 3667 observed with XMM-Newton . The main aim of 
the work is to study its metal distribution, that can re veal clues 
about the chemical enrich ment history of the cluster dSchindled 
2008; Borgani et al. 2008), and to infer its dynamical state by 
comparison with previous results and hydro dynamic simulations 
( Kapfe rer et al.ll2006l: IKapferer et al.ll2007|) . Secondly, we deter- 
mine the abundances of several elements and using the yields 
of Supernovae type la and II we try to understand the origin of 
metals in galaxy clusters. 

The paper is structured as follows: in Sect. 2 we give an overview 
of A3667; in Sect. 3 we present the data sets and data reduction 
techniques employed and we describe the X-ray image; in Sect. 
4 we present the metallicity map and relative comparison with 
simulations; in Sect 5 we present measurements of metal abun- 
dances and SN ratio determination, and in Sect 6 we discuss our 
results. A summary of our conclusion is given in Sect. 7. 
Throughout the paper we assume Ho=70 km s" 1 Mpc -1 , 
Q A =0.73 and Q M =0.27. At the nominal redshift of A3667 
(z=0.055), the luminosity distance is 245.7 Mpc and the angu- 
lar scale is 64 kpc per arcmin. 
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Fig. 1. X-ray contours (green) overlaid on the optical 
Digital Sky Survey image. The circles represent the main 
cluster (white full line) and the northwest subcluster 
(dashed line)', the ellipses represe nt the subgroup of 2 7 mem- 
bers (dashed yellow line) found by lOwers et al.l (2009) and the 
two radio relics (magenta lines). The image is 1° x 1°. 



2. Abell 3667 



A3667, a rich southern cluster at z=0.055 dSodre et al.lll992h . is 
famous for the presence of two extended radio relics symmetri- 
cally located in th e cluster periphery in th e direction of the elon- 
gated X-ray a xis jR ottgering et al ] 1 19971). It has a large veloc- 
ity di spersion dProust et alJl9 88. Sodr e et al.lll992llGirardi et al 



1998 ) and the 2D galaxy distribution is bimodal (jProust et al. 
1988L ISodre et al.lll992l ljohnston-Hollitt et alJl2008h , with the 



main component around the cD galaxy near the X-ray peak and 
the secondary component around the second bright est galaxy lo- 
cated in the northwest, -15' from the X-ray center. lOwers et alJ 
(2009), with a large number of confirmed galaxy members (550), 
found evidence for a new subgroup located in the southeast re- 
gion (see Fig. Q]). The bi modal structure is evident also in the 
weak lensing mass map (JoffreetaL|[2000) where it is possi- 
ble to see a significant mass concentration in the southeast of 
the cluster but not coincident with the substructure shown by 
lOwers et al] d2009b . ROSAT and ASCA observations have shown 
a distorted X-ray morpholo gy in the direction of the reported bi- 
moda l optical distribution (Kn opp et~aTlll996llMarkevitch et al] 
fl999h . XMM-Newton and Chandra observations have revealed 
an inhomogeneo us temperature structure dMazzotta et al. 2002, 



Briel et al .1120041) and evidence for a cold front (IVikhlinin et al] 
200 lb . All these features indicate that the cluster suffered a 



merger recently. Two different scenarios were suggested by 
Ower s et al. (2009). The first is a two-body merger between sim- 
ilar mass structures (the "main cluster" and the NW "subcluster" 
in Fig. [I]) taking place in the plane of the sky. In this scenario the 
NW subcluster would have already traversed the main cluster 
along a SE-NW direction, producing two outgoing shocks that 
would account for the double radio relics observed in A3 667. In 
this case, both the cold front and the subgroup in the southeast of 
Fig.[T]could have been associated with the northwest cluster, and 
then sloshed out or tidally stripped during the passage through 




Fig. 2. Total (MOS+pn) EPIC mosaic image of A3667 in the 
[0.5-8] keV. The image is non-background subtracted, corrected 
for vignetting and exposure and adaptively smoothed. Contours 
indicate the surface brightness spaced by a factor of V2 start- 
ing at 5.5 xlO -5 cnts s" 1 arcmin -2 level, which is 3.7 times the 
average background level. North is up and east is to the left, as 
they are in all other images. Coordinate grids are shown for the 
J2000. 



the core of the main cluster. The SE subgroup could also be a 
background or foreground structure. An alternative scenario in- 
volves a three-body merger between the main cluster and the 
NW and SE substructures along a NW-SE axis. The NW and 
SE radio relics would then be associated to the merger between 
the main cluster and the NW and SE subclusters respectively, 
and the cold front could be the remnant cold core of the SE sub- 
group. 

3. Observations and data reduction 

3.1. X-ray analysis 

Observation data files (ODFs) were retrieved from the XMM 
archive and reprocessed with the XMM-Newton Science 
Analysis System (SAS) v7.1.0. We used tasks emchain and 
epchain to generate calibrated event files from raw data. 
Throughout this analysis single pixel events for the pn data 
(PATTERN 0) are selected, while for the MOS data sets the 
PATTERNS 0-12 are used. In addition, for all cameras events 
next to CCD edges and next to bad pixels were excluded 
(FLAG==0). 

The data were cleaned for periods of high background due to the 
soft proton solar flares using a two stage filtering process. We 
first accumulated in 100 s bins the light curve in the [10-12] keV 
band for MOS and [12-14] keV for pn, where the emission is 
dominated by the particle-induced background, and exclude all 
the intervals of exposure time having a count-rate higher than 
a certain threshold value (the chosen threshold values are 0.20 
cps for MOS and 0.25 cps for pn). After filtering using the good 
time intervals from this screening, the event lists was then re- 
filtered in a seco nd pass as a safety check for possible flare s 
with soft spectra dNevalainen et al.ll2QQ5b iPradas & Kerpll2005h . 
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In this case light curves were made with 10 s bins in the full 
[0.3-10] keV band. The resulting exposure times after cleaning 
are 56.6 ks for MOS1, 56.1 ks for MOS2 and 45.7 ks for pn. 
To correct f or the vignetting effe ct, we used the photon weight- 
ing method (lArnaud et al.ll200l1) . The weight coefficients were 
computing by applying the SAS task evigweight to each event 
file. Point sources were detected using the task ewavelet in the 
energy band [0.3-10] keV and checked by eye on images gen- 
erated for each detector. We produced a list of selected point 
sources from all available detectors and the events in the corre- 
sponding regions were removed from the event lists. 
The background estimates were obtai ned using the dedicate d 
blank- sky event lists accumulated by Rea d & Ponmanl ( 2003). 
The blank- sky background events were selected by applying 
the same PATTERN selection, vignetting correction, flare re- 
jection criteria and point source removal used for the observa- 
tion events. In addition, we transformed the coordinates of the 
background file such that they were the same as for the associ- 
ated cluster data set. The background subtraction was performed 
using the double s ubtraction process described in full detail in 
lArnaud et al.l (120021) . It involves subtraction of the normalised 
blank field data, and subsequent subtraction of the cosmic X-ray 
background component estimated in the area of the field of view 
that does not show cluster emission. 

3.2. X-ray image 

As the X-ray morpholo gy can give intere sting qualitative (and 
quantitative, see e.g., iBuote & Tsailll996l) insights into the dy- 
namical status of a given cluster, the adaptively smoothed, ex- 
posure corrected (MOS+pn) count rate image in the [0.5-8] keV 
energy band is presented in Fig. [2 The smoothed images was 
obtained from the raw image corrected for the exposure map by 
running the task a smooth set to desired signal-to-noise ratio of 
20. Regions exposed with less than 10% of the total exposure 
were not considered. It is possible to see the elongated struc- 
ture of the cl uster and the co l d fron t to the southeast extensively 
discussed bv lVikhlinin etaD ([2001). The distortion in the north- 
west is introduced by the field of view (FOV) edge. In Fig. Q] we 
show the X-ray contours superposed on the optical image of the 
cluster. The X-ray peak lies at RA 20: 12.:27 and DEC -56:50: 1 1 
(J2000) and is near the central domin ant cluster galaxy located 
at RA 20:12:27 and DEC -56:49:36 (lOwers et alJl2007h - 

4. Metallicity map 

To obtain a metallicity measurement with a good accuracy re- 
quire a high statistic. Based on previous metallicity study we 
set a minimum count number (-3,000 source counts per region) 
necessary for proceeding with the spectral fit. The spectral re- 
gions for the map were selected using the following method. We 
first produced an image where each pixel is 500x500 EPIC phys- 
ical pixels corresponding to 25 // x25 // . So, from now on "1 pixel" 
is actually this "fat" 25"x25" "pixel". A square region with side 
length of -1100 arcsec and centred on the peak of the X-ray 
emission was defined to include only areas where the surface 
brightness of the source is high. This region was divided into 
11x11 square regions, each 100 arcsec 2 . The size of these square 
regions was then optimized by splitting it into horizontal or ver- 
tical segments through its center, while including at least 3,000 
source counts, summed over the all three EPIC cameras. Any 
region which did not contain 3,000 counts was ignored. For all 
the selected regions, spectra were extracted for source and back- 
ground in all three cameras. Finally, the spectra were re-binned 
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Fig. 3. Example of one EPIC pn (green) MOS1 (black) and 
MOS2 (red) spectrum used to determine the abundance for 
metallicity map with an error of 10-20%. 

with the grppha task, to reach at least 2 counts per en ergy bin. 
Spectra were analysed with XSPEC (Arnaudl Tl996l) version 
12.3.1. Since the spectra were re-binned, we have used stan- 
dard^ 2 minimization. We determined the errors with the XSPEC 
tasks error and step par. In order to model the emission from a 
single temperature we fit the spectra with the following model: 

MODEL = WABS (N H ) X MEKAL(T, Z, K). (1) 

WABS is the pho t oelectr ic absorption model by 
iMorrison & McCammonl (1 19831) and MeKaL model is the 
traditional plasma c ode dMewe et al .1 1 1 9 8 5L 1 1 9 8 6l: I Kaas tr al [19921: 
iLiedahl et al.lll995l) in which the temperature T, the metallicity 
Z and the normalization K are free parameters. The spectral 
fit was done leaving the hydrogen column density to free vary. 
We fit jointly MOS1, MOS2 and pn spectra, enforcing the 
same normalization value for MOS spectra and allowing the 
pn spectrum to have a separate normalization. In the spectral 
fitting we used the 0.5-8 and 0.5-7.5 keV energy range for MOS 
and pn spectra respectively. We excluded the energy above 
7.5 keV in the pn spectra because of the strong fluorescence 
lines of Ni, Cu & Zn. These lines, present in the background, 
are not well subtracted by the double background subtraction 
because they do not scale perfectly with the continuum of the 
particle-induced background. The redistribution and ancillary 
files (RMF and ARF) were created with the SAS tasks rmfgen 
and arfgen for each camera and each region that we analysed. 
The metal abundan c es are based on the solar values given by 
lAnders & Grevessel (Il989l) . We chose to use these abundance 
for easier comparison with previous work, although these values 
may not be accurate as shown by iGrevesse & Sauvall (1 19981) 
who obtained sign i ficantl y lower O and Fe abundance then 
lAnders & Grevessel (Il989h . In Fig. [3] we show an example of the 
spectra. 

The obtained metallicity and temperature maps are shown in 
Figs. [4] and [6] (w/?/?er panel). The resolution of the maps is the 
same (although it would be possible to obtain a good estimate of 
the temperature within smaller region) for a direct comparison 
between the two maps. Regions in white are those where the 
spectral signal to noise ratio was not sufficient to determine 
the metallicity. Both maps show a complex substructure as 
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Fig. 4. U pper panel. Metallicity map based on spectra from all 
three EPIC camera. The scale for the metallicity is in solar units. 
Lower panel. Number of bins with a certain metallicity. The ver- 
tical line in blue represents the mean value. 
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Fig. 5. Metallicity map for lcr upper (top) and lower (bottom) 
limits, displaced in the same color scale. The numbers indicate 
the four selected regions used for the abundances determination. 



expected for a merging cluster. In particular the metallicity 
distribution appears very inhomogeneous, while the temperature 
map shows a hot arc-like structure around the cold gas region. 
The highest peak of the metallicity is located in the southeast 
with respect to the X-ray center corresponding to the cold front 
region. A higher metallicity is also observed in two blobs to the 
NW. Between those clumps we note a region with a very low 
metallicity (below 0.2 in solar abundances). 
The histograms in Fig. [4] and [6] (tower panel) shows the distri- 
bution of the metallicity and temperature values. Concerning the 
metallicity we note that all the values range between 0.05 and 
0.7 solar abundances with a mean of -0.31 Z© that is in good 
agreement with the value of 0.29 obtained by fitting the spectra 
of the whole cluster within a radius of ~8 arcmin centered on 
the emission peak. 

Typical errors in our metallicity map are about 10-20%, al- 
though for a few pixels in the outskirts the error is higher. In 
Fig. [5] we show the lcr upper and lower limits for comparison. 
For the temperature maps all the errors are lower than 10%. 



5. Abundances and enrichment by supernova 
type la and II 

The chemical evolution of galaxies and consequently of the ICM 
is dominated by its main contributors SN la, SN II and plane- 
tary nebulae (the contribution of planetary nebulae is negligible 
for abundances of elements from O to Ni). We investigated the 
relative contribution of the supernova type Ia/II to the total en- 
richment on the intra-cluster medium. We assumed that the total 
number of atoms N; of the elements i is a linear combination of 
the number of atoms Yj produced per supernova type la (Yij a ) 
and type II (Y UI ): 



Ni = aY Ua +fiY iJ i, 



(2) 



where a and j3 are the numbers of supernova types la and II 
respectively. We used SNIa yields obtained from two different 
models adapted from llwamoto et all (1 19991) : the W7 model is a 
so-called slow deflagration mod el, while WDD2 is the favored 
model by llwamoto et al.l (1 19991) and is calculated using a de- 
layed detonation model. For nucleosynthesis products of SN II 
we adopted average yields of stars in a mass range from 10 M© to 
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Fig. 6. U pper panel: Temperature map based on spectra from all 
three EPIC cameras. Lower panel. Number of bins with a certain 
temperature. The vertical line in blue represents the mean value. 



50 M© calculated by Tsujimoto et a l. (1995) assuming a Salpeter 
initial mass function (IMF). 

From observations we obtained the number of nuclei per hydro- 
gen nucleus relative to the solar abundances: 



N H 



(3) 



where N# is the number of hydrogen atoms and fi the abundance 
obtained directly from XSPEC analysis. Combining equations [2] 
and [3] and considering iron as fixed element we can compute the 
ratio between supernova type la and II: 



a 



(Fe_ _ Ni\ 

\fFe fi), 



SNII 



(Nj_ _ Fe\ 
\ fi fFe tsNIa 



(4) 



where Fe in the equation is the yield value for the iron. We chose 
to fix iron because is the best determined element. 
We wanted to determine if the high metallicity region is associ- 
ated with a particular kind of SN type or with high number of 
SNe II as consequence of an intense e pisode of star format ion, 
due for instance to ram-pressure effect ( Kapferer et al. 2009|). To 
do this we selected four regions: first we extracted the spectra 
from a circular region, centered on the X-ray peak, with a radius 
of 3' (region 2 in lower panel of Fig. [5]) ; then we selected an 
annulus with inner and outer radius of 3' and 7' respectively and 




MOS 



Fig. 7. Upper panel: The hydrogen column density value ob- 
tained with a MeKaL model for the 4 considered regions us- 
ing the all three EPIC cameras ( MOS+pn). The red dashe d 
line represent the value obtained by lDickey & Lockmanl (Il990l) . 
Lower panel: Oxygen values obtained by fitting spectra from the 
MOS and pn detectors independently with the column density 
fixed to the value computed with a MeKaL model. 



we extracted a spectrum for the low metallicity region by select- 
ing a sector between PA = 1 80° and 90° (East to West) and for 
the high metallicity region selecting the sector between PA = 90° 
and 180° (respectively region 3 and 4); finally, we extracted the 
spectra from a circle with a radius of 7' (region 1) that includes 
all the three previous regions. 

5.1. Abundances determination 

In Table [T] we show the obtained abundances with their ±l<x 
errors for one parameter for the four selected regions. We fit- 
ted the data with the following procedure to avoid the degener- 
acy of the parameters: (1) we fitted the data with an absorbed 
MEKAL model in the 0.4-7 keV band to obtain temperature 
and nH (metallicity and normalization are considered free pa- 
rameters); (2) we fixed nH and temperature and use a VMEKAL 
model in the same energy band to determine the iron abundance 
(O, Si, S, Ar are left free, the other elements are fixed to the solar 
value); (3) we kept temperature and iron fixed to measure oxy- 
gen abundance in the 0.4-1.5 keV band; (4) we fix the values of 
temperature, iron and oxygen to estimate the silicon, sulfur and 
argon abundances in the 1.5-5 keV band. 

We fitted the element abundances in narrow bands around to 
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Table 1. Abundances obtained by fitting the four selected re- 
gions shown in Fig. [5j Temperatures are given in keV, N# in 
10 20 cm -2 and abundances are given with respect to solar. All 
errors quoted are at the 68% level for one interesting parameter 
(A* 2 =l). 



Par. 


1 


2 


3 


4 




N„ 


4.01+0.08 


3.95+0.09 


4.01+0.11 


4.08+0.12 




kT 


5.919+0.035 


6.157+0.053 


6.516+0.068 


5.050+0.049 


\_ 1.6 





0.34+0.05 


0.26+0.06 


0.30+0.07 


0.40+0.07 




Si 


0.43+0.07 


0.37+0.08 


0.35+0.07 


0.52+0.08 


1 


S 


0.18+0.10 


0.21+0.10 


0.12+0.10 


0.22+0.10 




Ar 


0.28+0.11 


0.27+0.18 


0.24+0.20 


0.36+0.20 




Fe 


0.306+0.009 


0.288+0.010 


0.275+0.014 


0.369+0.013 


0.5 


X\ed 


1.58 


1.19 


1.23 


1.23 





Table 2. Relative number of SN II contributing to the enrichment 
of the intra-cluster medium. We show the results for different re- 
gion extractions using the yields of SN la and SN II, with two 
physically different SN la yield models. In the second column 
we indicate the element used in combination with iron to deter- 
mine the percentage number of SN II. 



Model 



W7 



WDD2 



el. 



1 



O 
Si 
S 

Ar 



72 6 +3 - 7 
75.0^-6 

<55.4 



66.5^-1 

70.3^ 

<70.4 
84+15 

-84 



72.0^ 
69.8* 
<36.6 



71.9^ 
75.1^-3 

<48.8 
86 + g 



O 
Si 

s 

Ar 



75.5!^ 

<23.8 
75 +19 



70.0^ 75.0^ 



n+5.5 
7.4 

69.0 + _ l 2 \ 
<59.5 

IT 11 



+5.2 



74 ' 9 -5.1 



+7.0 
11 



7 Q+30 



74.5 
<2.1 



the corresponding emission lines, allowing the normalization to 
vary, in order to correct for small inaccuracies in the best deter- 
mination of the continuum in those narrow energy band. 
The Ne and Mg abundances could not be constrained because 
these lines are blended with the Fe-L complex at the EPIC spec- 
tral resolution. The aluminium line is blended with the much 
stronger silicon line and is not measurable. The Ni abundance 
determinations are driven almost entirely by the He-like and H- 
like K- shell lines at 7.77 and 8.10 keV which are both beyond 
the chosen spectral fitting band. 

In ¥\g.^\{upper panel) we show the hydrogen column density 
obtained in all the four analyzed regions. We note that is always 
lower than the Galactic valu e of 4.7 xlO 20 cm" 2 determi ned from 
the 21 cm radio observation iDickey & Lockmanl ([l990). We left 
free nH (not fixed to the Galactic value) because of the large 
discrepancy and also because the O abundance determination is 
sensitive to the presence of excess absorption and to the cross- 
calibration uncertainties between the spectral response of the 
two EPIC instruments in the soft energy band (below 1 keV). 
In the lower panel of Fig. [7] we show the oxygen abundance in 
the four different regions obtained by fitting the spectra from the 
MOS and pn detectors separately with N# fixed at the value ob- 
tained with a MEKAL model in the 0.4-7 keV band. The agree- 
ment between the two detectors is quite good; the discrepancy in 
the central 3' (region 2) can be due to a calibration problem at 
this particular position of the detector. 



Oxygen 
Argon 
Sulfur 
Silicon 



// / 




50 60 
TI (%) 

Fig. 8. The lines represent the abundance ratios as the SN ra- 
tio varying for the different elements. The two solid triangles 
present the observed O/Fe (red) and Si/Fe (black) ratio for the 
7' region while the boxes present the error associated with the 
observed abundance (Y axis) and the relative number of SN II 
associated with the lower and upper limit (X axis). 



6. Discussion 

6.1. Metallicity map 

There are a number of features to note analyzing the metal map. 
First, the metal distribution is very inhomogeneous with several 
maxima and complex meta l patterns as expected for a merging 
cluster dKapferer et al.ll20 06). The X-ray peak is located in a re- 
gion with a very low metallicity (0.15-0.20 Z , see Fig. [4]) with 
respect to the mean metallicity of the cluster (0.31 Z ). From 
the simulations is clear that the maximum of the m etallicity is 
not always in the cluster center (Kapf erer et al .1120061) . The rea- 
son is, that enriched gas, that might fall into the center, is mixed 
with a lot of other gas as in the centre the gas density is high. 
Therefore it hardly increases the metallicity there. If, however, 
a starburst happens in the outer parts of the cluster, the enriched 
gas mixed only with a small amount of other gas and therefore 
it can increase the metallicity there considerably. Hence it can 
happen that the maximum of the metallicity is temporarily not 
in the cluster centre. 

Althouh several blobs with high metallicity (-0.5 Z ) are present 
in the NW direction along the axis of the X-ray elongation we 
found that the most metal rich zone correspond t o the cold front 
region, as previously shown by lBriel et al.l (120041) . with a peak of 
0.7 solar abundance. This region probably shows the most inter- 
esting feature in the metal map and it can be used together with 
simulations to infer the dynamical state of the cluster. 

6.2. SN enrichment 

In Table [2] we show the obtained best fit values of the relative 
contribution of SN II with a confidence level of 68%. We see 
that both models are consistent with a scenario where the rela- 
tive number of supernovae type II contributing to the enrichment 
of the intra-cluster medium is 55-95% depending on the consid- 
ered elements and regions. The best agreement between the data 
of O, Si, Ar and Fe is obtained using the WDD2 model, although 
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the error bars are quite large due to the uncertainties both in the 
observations and in the theoretical yields. The relative number 
of SN II seems to be higher in the metallicity peak region (4), 
and lower for the regions 2 and 3. The lowest value is obtained 
in the center (region 2) and it is c onsistent with the i dea of an 
excess of SN la in the cD galaxies (iWerner et al.ll2008l) . 
The measured abundances of S are quite low and also consid- 
ering the lcr upper limit the percentage of SN II that we de- 
rive is not consistent with the results given by the other ele- 
ments. We note that the relative low value of 0.18 obtained for 
the abundances of S in a radiu s of 7 ' is in good agreement with 
the value of 0.20 computed by iBriel et al.l (120041) for the whole 
cluster. This value agrees also with the result of the sulfur abun- 
dance obtained with ASCA d ata for a sample of clust ers with the 
same temperature as A3 667 teaumgar tner et al.ll2005l) confirm- 
ing both the prevalence of SN II in the enrichment and a reduced 
S yield in the SN II model. 

We note that for all the elements the lower relative error is larger 
than the upper one in the relative SN II determination. To explain 
the reason for it, we combine the equations [2] and [3 and we find 



f Fe Fen + FeA' 



(5) 



that gives the theoretical abundance ratio ft I fFe for varying SN 
type Ia/II ratios. We show in Fig. [8] the ft I fFe obtained using the 
yields of the WDD2 model for supernovae type la. It is clear 
that not all the values of fl fFe are allowed: for example the 
Ar/Fe abundance ratio must be in the range between 0.7 and 1 .5. 
In Fig. [8] we see that all the curves grow very slowly for a low 
percentage of SN II and they become very steep when that per- 
centage increases. This implies that the determination of the SN 
fraction becomes quite inaccurate when the relative number of 
SN II is low, in particular for S/Fe and Ar/Fe that shows a very 
flat curve when the number of SN II is lower than 60% of the 
total supernovae. If we consider symmetric error bars for fl fFe 
it is clear why for the lower relative error in the SN II determina- 
tion is always larger than the upper one. This does not influence 
too much the determination of the SN relative number only if the 
curve is very steep as for example the case of the O/Fe ratio. 
To obtain the best agreement between the observed abundances 
of O, Si, Ar and Fe for the 7' region we estimated that 65-80% 
of SN II are necessary. Other authors tried to determine the ratio 
of SN la to SN II events in relaxed galaxy clusters by aiming 
for a best fit to a n ove rall solar abundance pattern from O to 
Ni. I Werner et al.l (120061) found that the number contribution of 
SN II with respect t o the total number of supernovae is -75%, 
Ide Plaa et al] (l2006|) constrained th is number to the range 50- 
75% while ISimionescu et al.1 d2QQ9h estimated a 30-40% contri- 
bution by SN la compared to type II. These results suggest that 
so far th e enrichment of the ICM is mainly due to SN II. We 
note that lPipino et al] (120021) using different chemical evolution 
models for galaxies, showed that SN II dominate the chemical 
enrichment inside the galaxies, while la supernovae play a pre- 
dominant role in the ICM, that is not in agreement with the ob- 
servational results. 

When we interpret the supernovae ratio we have to take into ac- 
count that the abundances ratio does not only depend on stel- 
lar yields and IM F but also on the timescales of production of 
various elements ( Mat teucci & C hiappini 2005). The abundance 
ratios will tend to the ratios of their yield per stellar generation 
only if the global metal production is considered (metals in stars, 
gas inside and outside the galaxies), but it fails if only the metals 
in the individual component are taken into account (e.g. the gas 



of ICM). Thus, the supernovae estimation listed above should be 
interpreted as the number of supernovae that would be needed 
to reproduce the same abundances observed in the ICM, and not 
the number of supernovae during the history of the cluster. 
It is interesting to compa re the est i mated values with the ones 
obtained for the galaxies. iLeamanl (120081) shows the results of 
the Lick Observatory Supernova Search (LOSS) and in partic- 
ular he focuses on the determination of the supernovae in the 
local universe. He found that about 62% of the SNe observed 
in gala xies are SN II. In order to reproduce the observed abun- 
dances, |Ts^imotoetaL| (1 19951) determined the percentage of SN 
II for our Galaxy to be -87%, while it is -77% and -83% for 
the Large and Small Magellanic Clouds respectively. The rela- 
tive contribution to the enrichment of ICM by SN II in A3 667 is 
between these values. 



6.3. Comparison to simulations 

Simulations of the enrichment of the ICM are a powerful 
tool to investigate the strengths of different enrichment pro- 
cesses and their spatial and temporal behavior. To study the 
origin of the inhomogeneities found in the metal distribu- 
tions in many galaxy clusters, several simulations were per- 
formed. In the numerical setup we used different code mod- 
ules to calculate the main components of a galaxy clus- 
ter in the framework of a standard ACDM cosmology. The 
non-baryonic da rk matter (DM ) component is calculated us- 
ing GADGET2 ( Springel 2005) with constr ained random field 
initial conditions (Hoffman & Ribak |jl99ll) . implemented by 
Ivan de Weygaert & Bertschingerl (ll996l) . For the treatment of the 
ICM we use comoving Eulerian hydrodynamic wit h a shock cap- 
turing schemes ( PPM. IColella & Woodwardlfl984l) . with a fixed 
mes h refinemen t ( Ruffert 1992) on four levels and radiative cool- 
ing dSutherland & Dopitalll993l) . The properties of the galaxies 
are calculated by an improved version (Ivan Kampen et al.ll2005l) 
of the galaxy formation code of Ivan Kampen et al] (1 19991) which 
is a semi-analytic model in the sense that the merging history 
of galaxy halo is taken directly from the cosmological N-body 
simulation. With this setup we investigated two different enrich- 
ment processes, namely su pernova driven galac tic winds and 
ram-pressure stripping (see KapfeTerl let al. 2009| for details re- 
garding the galactic winds model and lDomainko et ai]|2006l for 
the ram-pressure stripping model). 

In Fig. [TO] the evolution of the metallicity in a model cluster is 
shown. The isolines correspond to the X-ray surface brightness 
of the model cluster. In the simulation a region of enriched ma- 
terial falls towards the cluster center. The high metallicity is in 
this infall region (see at the top of panel (a) in Fig.[l0|) and it is 
caused by four starbursts (with an outflow rate of more than 25 
solar masses each) that happened before z=0.42. Approaching 
the cluster center the material gets mixed with the lower metal- 
licity ICM along the trajectory (see panel (b) in Fig.[l0|). From 
the turnaround point (see panel (c) Fig. [TO]) the material falls 
smoothly to the center and gets mixed by the ambient ICM with 
lower metallicity (see panel (d) Fig.[l0|). As the material ejected 
by galactic winds and starbursts contains more SNII products 
the feature in the simulation (panel (d) in Fig. [TOj) corresponds 
nicely to the off-center metal concentration found in A3 667 (see 
Fig. [4]) where the relative contribution of the SN II is higher (re- 
gion 4). At a redshift of z=0 the enriched gas originating from 
starbursts has already mixed with the ICM, leading to the metal 
map shown in panel (d) in Fig.fTOl 

Metal blobs originating either from galactic winds or ram- 
pressure stripping are a common feature in metal enrichment 
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Fig. 9. Metallicity maps with the galaxies isodensity (left) and X-ray (right) contours superimposed. The scale of the metallicity is 
in solar units. An elongated substructure, both in optical and X-ray contours is shown. 



simulations. Typically they move along the trajectory of the un- 
derlying galaxies and as they start to feel the pressure of the 
surrounding gas they lag behind the originating galaxies and 
mix with the surrounding gas. The more inhomogeneities are 
present in the ICM the more recent the enrichment processes 
took place. Therefore the inhomogeneities found in the metal 
maps in galaxy clusters are indicators for the merging frequency 
of substructures with the cluster. Typically the inhomogeneities 
in the metal maps vanish over timescales of several 100 Myrs to 
several Gyrs depending to the mass present in the metal feature. 
In the example in Fig.[l0|the metallicity blob survives nearly 3 
Gyrs. 

Based on the size of the high metallicity region (~3' of radius), 
we expect that the metal feature in Fig.|4]is either a consequence 
of a recent merger or of an older merger but which involve a 
larger mass. In the latter case the metal feature would have sur- 
vived for a long time. 



6.4. Dynamical state 



We p roduced a multi-scale galaxy density map (see lFerrari et al.l 
2005 for more details) using t he 550 spectroscopic ally confirmed 
cluster members obtained by lOwers et al.l (l2Q09h . In Fig. [9] we 
plotted for comparison the metal map with the galaxy isoden- 
sity overlaid (left) and X-ray (right) contours. As for the X- 
ray surface brightness the projected galaxy density map shows 
an elongation in the direction of the two radio relics. The com- 
parison of the ICM metallicity distribution and the position of 
the sub clusters of A3 667 c an give hints on the c omplex merg- 
ing scenario of this cluster dKapferer et al.ll2006h . We detected 
a metal peak between the main cl uster and the SE subgroup. 
According to iKapferer et al.l (120061) we expect high metallicity 
between subclusters in a post-merger phas e. Thus, this config- 
uration supports the scenario suggested by lOwers et all d2009h 
in which the SE subgroup has traveled from the NW and passed 
through the main cluster where the ram pressure stripped off the 
enriched and cooler gas. 



However, an elongation towards the high metallicity peak visi- 
ble both in optical and X-ray images (see the contours in Fig. 
[9]) could suggest a more complex dynamics in the cluster center. 
The abundances of the measured elements are higher in region 4 
with respect to the regions 2 and 3. This result can be explained 
if a group of galaxies, located in the elongated substructure and 
containing both SN la and SN II products, falls into a cluster 
moving from southeast to northwest. The inter- stellar medium is 
thus stripped off by ram-pressure stripping and leads to a more 
peaked abundance distribution. On the other hand, in region 4 
the relative number of SN II seems to be higher, with respect 
to the other two considered regions (region 2 and 3), suggesting 
that the metallicity peak in region 4 is mainly due to galactic 
winds as obtained in the simulations. Due to the large error bars 
in the SN determination the latter result has to be confirmed with 
a deeper observation. 



Inter estingly a gap in metallicity has been detected by lBriel et al.l 
2004 in between the main cluster and the NW sub-cluster (i.e. a 
regi on not covered by o ur metallicity map). Based on the results 
of dKapferer et al.ll20 06) this would imply that these two struc- 
tures are in a pre-merger phase. All these results could suggest 
that A3 667 is a cluster forming through multiple merging events 
along a common NW-SE axis. 

7. Summary 

We analyzed a 64 ks XMM -Newton exposure of the merging 
cluster of galaxies A3 667. We obtained a detailed 2D metallicity 
map. From this we can conclude that: 

- the distribution of metals is clearly non- spherical. It looks 
very inhomogeneous with several maxima separated by very 
low metallicity regions; 

- the highest metallicity peak is located on the southeast with 
respect to the X-ray center and it corresponds to the region 
with the lowest temperature. 

We also measured the abundances for oxygen, silicon, sulfur, ar- 
gon and iron in 4 different regions of the cluster and determined 
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the number ratio of supernovae type la and type II. From these 
data we conclude that: 

- using the elements abundance of Fe, O and Si we found that 
the relative number of supernovae type II necessary to re- 
produce the observed abundances in A3 667 ranges between 
65-80%; 

- the delayed detonation model WDD2 seems to reproduce 
the observed data better compared to the slow deflagration 
model; 

- the supernovae number estimation from the abundances of 
sulfur is not in agreement with the estimate obtained using 
the other elements confirming a reduced S yield in the SN II 
model. 

Finally, we discussed the dynamical state of the cluster by com- 
paring the ICM metal and galaxy density maps to our simula- 
tions. In agreement with the scenario proposed by lOwers et al.l 
(2009), we conclude that the SE subgroup moved from the NW 
and passed through the main cluster, where the ram pressure 
stripped off the enriched and cooler gas as seen in the metal- 
licity and temperature maps. The highest metallicity region, that 
shows a higher contribution of SN II, could be partly related to 
an enrichment by galactig winds due to star formation possi- 
bly triggered by an inf ailing group. In addition, two metal rich 
blobs in the NW of the main cluster could partly result from 
inhomogeneities not completed dispersed after an old merger, 
which is possibly responsible for the formation of the two radio 
relics. Based on the comparison with previous X-ray data (e.g. 
iBriel et al.ll2004 we conclude that A3 667 has a complex dynam- 
ical history and it is possibly evolving by accreting sub-clusters 
along a main NW-SE axis. 
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L. Lovisari et al.: Metallicity map of the galaxy cluster A3667 




Fig. 10. Simulated X-ray weighted metal maps. X-ray surface 
brightness contours are overlaid. A model cluster showing sim- 
ilar features in the metallicity distribution as A3 667 as been 
selected. The five maps correspond to different redshifts: (a) 



